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ABSTRACT 

We present a map of the Cosmic Microwave Background (CMB) anisotropics 
induced by the late Integrated Sachs Wolfe effect. The map is constructed by combining 
the information of the WMAP 7-yr CMB data and the NRAO VLA Sky Survey 
(NVSS) through a linear filter. This combination improves the quality of the map 
that would be obtained using information only from the Large Scale Structure data. In 
order to apply the filter, a given cosmological model needs to be assumed. In particular, 
we consider the standard ACDM model. As a test of consistency, we show that the 
reconstructed map is in agreemet with the assumed model, which is also favoured 
against a scenario where no correlation between the CMB and NVSS catalogue is 
considered. 

Key words: cosmology: observations - cosmology: large scale structure - methods: 
data analysis - methods: statistical 



1 INTRODUCTION 

Recent observations, as cosmic microwave background 
(CMB), supernovae type la (SNIa) or baryon accoustic os- 
cillations (BAOs) agree in establishing a current accelerated 



and the galaxy number density field. Posterior works (e.g., 



expansion of the Universe (see Weinberg et al. 2012 for a 
recent review), which, despite some less popular interpreta- 
tions, is believed to be caused by the presense of some dark 
energy (see |Peebles fe Ratra|20"03| for a review) . The nature 
of dark energy is one of the most puzzling issues in mod- 
ern cosmology. The actual characteristics of this fluid is still 
unclear, although, up to date, a good agreement is found 
between the observations and the predictions derived from 
the presence of a cosmological constant with an equation of 
state p — —p. 

One of the classical probes of dark energy is given by 
a non-null contribution to the CMB anisotropics from the 



late integrated Sachs- Wolf effect (ISW |Sachs fe Wolfe|1967) , 
under the assumption of a spatialy flat universe. The (lin- 
ear) ISW fluctuations are higher at very large angular scales, 
but, in any case, much smaller than the primary CMB fluc- 
tuations. In a seminal work, [Crittenden &; Tur ok ( 1996) pro- 
posed the cross-correlation of the CMB fluctuations and the 
dark matter distribution (typically traced by galaxy cata- 
logues) as a possible approach to detect the ISW. Soon after 



the release of the WMAP data, Boughn & Crittenden (2004) 



Fos alba et al.|2003||Vielva et aI.|2006||Pietrobon et al.| 2006 

Ho et al. I20081 IGiannantonio et al. 2008; McEwen et al. 2008 



Dupe et al.||2011| |Schiavon et al.||2012| ) have confirmed the 
detection of the ISW effect by exploring several galaxy cat- 
alogues and cross-correlation techniques. Average detection 
is found at ~ 3a. 

Besides these works focused on the statistical detection 
of the ISW, more recently there have been attempts to re- 
cover the actual ISW fluctuations on the sky. In principle, 
an optimal ISW map could be derived from a 3D gravi- 
tational potential. This has been explored from very large 
simulations (e.g., Cai et al.|2010 ). The complexity to recover 
optimally the potential from surveys of galaxy with redshift 



information is challenging (e.g., Kit aura et al.||2009| |Jasche| 
et al.|2010| ), but very promising from the ISW studies point 
of view (e.g., |Frommert et al.||2008| ). 



Other approaches can be followed from surveys where 
the redshift information is poor, or known only statisti- 
cally. |Barreiro et al. (2008) proposed to use jointly maps 
of CMB anisotropics and of the galaxy number density field 
to recover the ISW signal on the sky. Alternative works mak- 
ing use only of galaxy catalogue maps have been proposed 



reported the first detection of the ISW efefct via the CMB 



2MASS. 



afterwards: Granett et al. (2009) on LRGs from SDSS-DR6, 



Francis fc Peacock (2010) and Dupe et al. (2011) on 
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This paper presents an application of the approach de- 



scribed in |Barreiro et al.| (|2008|) to W MAP (Jarosik etltL] 



2Q11D and NVSS ( | Condon et al.|1998| ). The outline of the ar- 



ticle is as follows. The methodology is reviewed in Section [2] 
A description of the data used, as well as the fiducial theoret- 
ical model is presented in Section [3] In Section [4] we present 
the results. Finally, conclusions are given in Section [5] 



2 METHODOLOGY 

In order to reconstruct the ISW map, we have used the linear 



covariance-based filter presented in Barreiro et al. (2008) 



We give here the outline of the method. 

Since the filter is implemented in harmonic space, for 
simplicity, we will assume that the considered data sets are 
full-sky. Let us denote se m and ggm to the harmonic coeffi- 
cients of the ISW map and the large-scale structure (LSS) 
survey respectively. The covariance matrix C{£) of the sig- 
nals at each multipole £ is given by 



C(£) 



CI 

c 3 / 



Cf 

ct 



(1) 



where and C| correspond to the auto spectra of the 
galaxies and ISW maps, respectively, while C s g 9 is the cross 
power between both signals. To construct the filter, we will 
make use of the Cholesky decomposition of the covariance 
matrix, which satisfies C(£) = L(£)L T (£), where L{£) is a 
lower triangular matrix. It can be trivially shown that the el- 
ements of the Cholesky matrix relate to the elements of C (£) 



as Ln 



'C 9 t , L12 



V\c(i)\/cf, 

where |C(£)| is the determinant of the covariance matrix at 
each £ mode. 

The estimated ISW map S£ m at each harmonic mode is 
given by (see |Barreiro et aLp 008 for details) 
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(2) 



where dim are the harmonic coefficients of the CMB map 
and Ci is the power spectrum of the CMB signal without 
including the ISW. Therefore, to reconstruct the ISW map, 
we need to assume an underlying cosmological model that 
determines the auto and cross power spectra present in the 
previous equation. 

It is interesting to note that the final reconstructed map 
has two contributions: the first term in the previous equa- 
tion is given by a filtered version of the galaxies map while 
the second expression is a Wiener filter (WF, |Wiener|1949| ) 
of a modified CMB data map. This modified data is sim- 
ply constructed as the original CMB map minus the filtered 
survey. In the case that there is not correlation between the 
CMB and the LSS survey, the filter simply defaults to the 
WF of the CMB map: since there is not correlation between 
both signals, the galaxies map does not contribute to the 
final ISW reconstruction. If the information provided by the 
CMB map were not considered, the estimated ISW would 
be given just by the filtered galaxies map. 

It can be easily shown that the expected value of the 
power spectrum of the estimated ISW is given by 



(c s /f(\c(£)\ + c?c?) + \c(e)\ 2 

C°(\C(t)\+CfC?) 



It is well known that the power spectrum of the WF recon- 
struction is biased towards values lower than the true signal, 
with the bias depending on the signal-to-noise ratio of the 
data. Since our signal is partially reconstructed using this 
filter, it will also be biased. The larger the cross-correlation 
between CMB and the considered galaxies catalogue, the 
smaller the bias, since the WF part will contribute relatively 
less than the filtered survey term (see Barreiro et al.| |2008 
for details). 

It is also straightforward to show that the expected 
cross-correlation between the recovered signal and the galax- 
ies catalogue is equal to that of the assumed model. How- 
ever, this is only correct if the assumed cross and auto power 
spectra reflects the underlying statistical properties between 
the ISW and the LSS survey. In practice, the difference be- 
tween the expected values of the previous quantities for two 
(reasonably) different models will in general be small and, 
given the weakness of the signal and the statistical uncer- 
tainties, it may be difficult to discriminate between them. In 
any case, as we will see in Section [2] the comparison between 
the expected and estimated values of the cross-correlation is 
an interesting consistency check. 

The previous description assumes that full-sky data are 
available. However, in practice, a mask will be needed to ex- 
clude those regions, in both CMB and galaxies survey maps, 
that have not been observed or are too contaminated to 
be included in the analysis. In Ba rreiro et al.| ([2008 ) was 
shown that the method was robust against the presence of 
a mask and that the quality of the reconstruction was not 
significantly affected. Therefore, to deal with this problem, 
we will simply substitute in the previous equations the har- 
monic coefficients and power spectra by those obtained after 
masking the data with the considering mask. In particular, 
the masked version of the fiducial model for the power spec- 
tra will be obtained ala MASTER ( |Hivon et aLl|2002| ). In 
addition, we will make use of an apodised version of the 
considered mask in order to reduce the correlations between 
harmonic modes that are introduced on an incomplete sky. 



3 DATA DESCRIPTION 

In order to reconstruct the ISW map using the previous 
methodology, we need both a CMB map and a LSS cata- 
logue, as well as a cosmological fiducial model. For the latter, 
we have assumed a ACMD model that best fits 7-yr WMAP 
data, BAO and Ho measurements ( Komatsu et al.||2011 ). 

For the CMB, we h ave made use of the 7-yr WMAP 
data ( | Jarosik et al.|[20TT| publicaly available at the Legacy 
Archive for Microwave Background Data Analysis Q In par- 
ticular, we have constructed a CMB map as a noise- weighted 
combination of the foreground reduced V- and W-band 
maps. The map has then been downgraded to a HEALPix 
resolution of N S id e = 32. To reduce the galactic foreground 
contamination, the KQ85 mask provided by the WMAP 
team has been used. However, since we have to downgrade 
the mask to a lower resolution, we have not included in the 
mask the holes due to point sources since their contribution 
is expected to be negligible at the considered scales. The 



(3) 
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(2011 ), where the bias is just defined by the minimum halo 





Figure 1. Top: WMAP 7-yr data (in Kelvin thermodynamic tem- 
perature). Middle: NVSS galaxy density number field. Bottom: 
apodised mask used for the analysis. 



WMAP data, with the considered mask applied, is given in 
the top panel of Fig. [I] 

Regarding the LSS catalogue, we have used he NRAO 
VLA Sky Survey (NVSS, [Condon et al.||l998) . NVSS cov- 
ers the north hemisphere of the sky and part of the south 
until S = -40°. The NVSS catalogue has nearly 2 x 10 6 
discrete sources with fluxes above 2.5 mJy. The catalogue 
was made using two different configurations (D and DnC) 
of the Very Large Array (VLA). This introduced a decli- 
nation dependence in the number of sources, specially for 
the faintest ones. Sources above 5.0 mJy have been selected 
to build a catalogue not affected by this declination effect. 
The NVSS mask is defined by the unobserved sky, plus the 
sources within 7 degrees from the galactic plane, as well as 
some regions associated to very near objects. The complete 
mask associated to NVSS leaves 73% of the sky. The map is 
shown in the middle panel of Fig. [I] 

The ISW-NVSS angular cross-power spectrum is cal- 
culated using the galaxy redshift distribution proposed by 



mass associated to NVSS radio-galaxies. The value of this 
mass is determined by fitting the NVSS angular power spec- 
trum (see Ma rcos-Caballero et al.|2012 , in preparation, for 
details) . 

Hernandez-Monteagudo 2010), 



As it is well known (e.g. 



de Zotti et al.| (|2010|: a p olynomial fit to the CENSORS 
data ( Brookes et al.||2008 ). For the galaxy bias is assumed 
a redshift dependence as the one proposed by |Xia et aT] 



there is a clear discrepancy between the meassured NVSS 
auto power spectrum and the theoretical model as desbribed 
above. In particular, an excess of power at low multipoles is 
observed. This deviation does not seem to be related to the 
declination effect previously mentioned, since it is observed 
even for the angular power spectrum estimated from very 
bright sources. A possible explanation in terms of a non- 
linear evolution of the galaxy bias (caused by a primordial 
non-Guassiniaty of the primordial perturbations) has been 
considered by Xia et al. (2011), although it would imply a 
value of the non-linear coupling parameter /nl larger than 
the one obtained from other observables (as the CMB bi- 
spectrum estimated from WMAP; see, for instance, Curto 
et al.|2Qj"T| . Currently, the nature of this discrepancy is not 
clear. Since the filter described in the previous section (see 
equation [2]) is very sensitive to the value of the auto-power 
spectrum of the galaxy survey, we have used as fiducial Cf 
a smooth fit to the observed NVSS spectrum. 

Finally, we have combined the WMAP and NVSS 
masks, leaving a useful area of 66%. The mask is then 
apodised using a cosine function, in order to reduce the cor- 
relations among the harmonic modes of the incomplete sky. 
This apodised mask (see bottom panel of Fig. [I]) is applied 
to the CMB and NVSS data prior to the reconstruction of 
the ISW map. In addition, the monopole and dipole outside 
the combined mask are subtracted prior to the construction 
of the ISW map. 



4 ISW RECONSTRUCTION 

Following the described methodology, we have produced an 
ISW map from the WMAP and NVSS data. The map is 
given in Fig. [2] and has been constructed at a HEALPix 
resolution of N S ide = 32 and up to a maximum multipole 
I = 95. 

As a consistency test, we have obtained the cross cor- 
relation between the reconstructed ISW map and the NVSS 
data and compared it with the fiducial model. The former 
has been simply obtained as the pseudo cross spectrum from 
the masked maps and then corrected ala MASTER. In addi- 
tion, we have also repeated the same procedure on two sets 
of 10000 simulations. For the first set, we have constructed 
correlated CMB and NVSS-like simulations (generated as 
explained in Barreiro et al.|2008| ) , assuming the same power 
spectra as for reconstructing the ISW map. For the second 
set, the autospectra are the same but zero correlation is as- 
sumed between the CMB and NVSS simulations. In both 
cases, we have included the corresponding Poissonian noise 
in the simulations of the galaxies catalogue. The results for 
this test are given in Fig. [3] 

As expected, the average value (red triangles) for the 
correlated simulations agrees well with the considered model 
(black line) . More interestingly, the data (blue asterisks) are 
also compatible with these simulations, which provides a 
good consistency check for our procedure and indicates that 
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Figure 2. Reconstructed ISW map. Units are Kelvin (thermodynamic temperature). 



the assumed cross and auto spectra are reasonably close to 
the true statisitcal properties of the sky. 

Conversely, the average obtained from the uncorrelated 
simulations (green squares) are biased with respect to the 
assumed cross spectrum. This is expected, since the assumed 
model is not the same as the underlying true power spec- 
trum. Unfortunately, since the signal is very weak, the dif- 
ference between both models is also small and it is difficult 
to distinguish between both hypotheses in the basis of this 
test. 

To try to quantify which hypothesis is favoured by the 
data, we have calculated the goodnes-of-fit difference: 



Ax 2 
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2 

X^ 



= Xc 
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i = {c, u}, and c sg correspond to the vector constructed 
with the values of the binned cross power spectrum. M c 
is the covariance matrix of the cross-spectrum obtained 
from 10000 correlated simulations and M u the same ma- 
trix for the uncorrelated simulations set. We have calculated 
the distribution of this quantity from independent sets of 
10000 correlated and uncorrelated simulations and compare 
it to the value found for the data. In particular, we find 
Ax 2 (data) = —2.4. For the correlated CMB and galaxies 
survey simulations, we find that 28.9 per cent of the sim- 
ulations have a value larger than the one obtained for the 
data. This percentage is reduced to 6.7 per cent, when we 
compare the data with the distribution obtained for uncor- 
related simulations. Therefore, although the data are in fact 
compatible with the two cases considered, it favours the hy- 
pothesis of having an underlying correlation. This further 
confirms the consistency of our results. 



Although, given the weakness of the considered cross- 
correlation, the power of these tests is somewhat limited, 
the fact that the data shows consistency with the assumed 
model should not be taken for granted, since other assump- 
tions may also affect the results. For instance, if we use the 
theoretical model for the auto spectrum of the NVSS map, 
instead ot a smooth fit to the measured data (as explained 
in section |3| , we find a clear departure of the cross corre- 
lation between the reconstructed ISW map and the NVSS 
data with respect to the expected value. 

Using simulations, we have also studied which is the im- 
provement from combining both the CMB and NVSS data, 
with respect to use only the galaxies catalogue (as for in- 
stance in Du pe et al.||201l| whose proposed reconstruction 
corresponds to the term given by the filtered survey of equa- 
tion]^. In particular, we find an improvement of 15 per cent 
in the error of the ISW reconstruction when the CMB data 
is included. In any case, the relative contribution to the re- 
construction of the CMB and LSS data depend significantly 
on the cross correlation between them, so it will differ for 
other galaxies surveys. 

In a related matter, we can also calculate the expected 
relative contribution of the auto and cross power of the CMB 
and the galaxy number density maps to the angular power 
spectrum of the reconstructed ISW map (using equation 
This is shown in Fig. [4] for the fiducial model as well as the 
actual contribution of each angular power spectrum for the 
considered data. As can be seen, the contribution from the 
power spectrum of the CMB C|* is dominant at low mul- 
tipoles (up to t ~ 10), while C 99 gives the main weight at 
higher multipoles. The contribution from the cross power 
between the CMB and NVSS is subdominant at all multi- 
poles. The results obtained from the data agrees quite well 
with the expected value. However, it is interesting to point 
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Figure 3. The binned cross-correlation between the recon- 
structed ISW map and the NVSS data is given (blue asterisks). 
For comparison, the fiducial model (black line) and the average 
and dispersion from two sets of 10000 simulations are also shown. 
The red triangles and their error bars correspond to CMB and 
NVSS-like simulations which are correlated, while for the second 
set of simulations (green squares are error bars), no correlation is 
introduced. The width of the bins is = 9, except fot the last 
bin that is A^ = 13. For a better visualization, the results for the 
simulations have been shifted with respect to the central value of 
the bin. 
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Figure 4. The expected relative contribution of the CMB and 
NVSS auto and cross spectra to the power spectrum of the ISW 
reconstruction is shown (solid lines) as well as the same quantities 
for the data. The results for the data have been binned to allow for 
a clearer comparison. The error bars correspond to the dispersion 
of the values within the considered bin. 



data. In particular, we have used the 7-yr WMAP data and 
the NVSS galaxies catalogue. The joint combination of both 
data sets improve in a 15 per cent the error of the ISW re- 
construction in comparison to the case when only the NVSS 
data are used. We have performed a consistency test, show- 
ing a good agreement between the cross correlation inferred 
from the reconstructed ISW map and the assumed fiducial 
model. In particular, the data favours a ACDM model with 
respect to a scenario with null correlation between the CMB 
and the NVSS data. The relative contribution to the angular 
power spectrum of the ISW reconstructed map is dominated 
by the CMB fluctuations up to £ « 10 and by the density 
number galaxies field at larger multipoles. 

The presented methodology works in harmonic space, 
which implies the use of surveys with large sky coverage, in 
order to avoid the problematics introduced by large masks. 
However, this technique can be easily extended to work di- 
rectly in the pixel space, which would make straightforward 
to deal with a mask and, in particular, LSS surveys with 
smaller sky coverage could also be used to reconstruct the 
ISW signal ( |Bonavera et al.||2012[ in preparation). In addi- 
tion, several catalogues can be combined at the same time, 
provided the covariance matrix between the surveys and the 
CMB data are known. 

Finally, let us remark that the application of the ap- 
proach described in this paper to future surveys (as EU- 
CLID |Refregier et~aLl|20Tol or J-PAS |Benitez et aLl|2009l 
with very large sky coverage and very accurate redshift esti- 
mation) could provide maps of the ISW anisotropics caused 
by the large-scale structure at different redshift shells. This 
will provide a tomographic view of the ISW fluctuations. 
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out that the relative importance of C|* with respect to Cf 9 
for the data is lower than the expectation value. This can 
be easily explained by the fact that the amplitude of the 
low multipoles of the WMAP data is lower than that of the 
fiducial model. 



5 CONCLUSIONS 

We have presented a reconstructed map of the ISW effect, 
obtained by combining through a linear filter CMB and LSS 
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